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We discuss the behavior of two point functions for 
Polyakov loops in a SU(3) gauge theory about the critical 
temperature, T c . From a Z(3) model, in mean field theory 
we obtain a prediction for the ratio of masses at T c , extracted 
from correlation functions for the imaginary and real parts of 
the Polyakov loop. This ratio is mi/m r — 3 if the potential 
only includes terms up to quartic order in the Polyakov loop; 
its value changes as pentic and hexatic interctions become im- 
portant. The Polyakov Loop Model then predicts how mi/m r 
changes above T c . 



Consider the behavior of an SU (3) gauge theory, with- 
out dynamical quarks, at nonzero temperature. The 
usual quantity measured by the Lattice is the pressure. 
While this is all one needs for thermodynamics in equilib- 
rium, near equilibrium how correlation functions behave 
is also of importance. 

In this note we discuss a basic, gauge-invariant cor- 
relation function, the two-point function of Polyakov 
loops. We work with a Z(3) model |],|; near the criti- 
cal temperature, this is equivalent to the Polyakov Loop 
Model |-§. 

At a nonzero temperature T, the thermal Wilson line 

is: 



L(x) = "Pexp ig 



l/T 



Aq (x, t) dr 



(1) 



with Aq in the fundamental representation. This trans- 
forms as an adjoint field under the local SU(3)/Z(3) 
gauge symmetry, and as a field with charge one under 
the global Z(3) symmetry. To obtain a gauge invariant 
operator, the simplest thing to do is to take the trace, 
forming the Polyakov loop, 



t = itr(L) 



(2) 



This transforms under the global Z(3) symmetry as a 
field with charge one, 



exp | — — ] £ 



(3) 



The expectation value of I is only nonzero above T c , 
which is the temperature for the deconfining phase tran- 
sition. 



For three (or more) colors, I is a complex valued field. 
Implicitly, we assume that the expectation value of £ is 
real. In this vacuum, we can then ask for the correlation 
functions of the real and imaginary parts of the Polyakov 
loop, 



Re^ 



(4) 



which do not mix. As the theory is globally Z(3) invari- 
ant, we can define equivalent masses in either of the other 
two degenerate vacua by an appropriate overall Z(3) ro- 
tation. 

We need a standard feature of propagators in coor- 
dinate space. In d space-time dimensions, at distances 
large compared to l/T, only the d — 1 spatial dimensions 
matter. For a field of mass m, at distances x 3> 1/m, of 
course the propagator falls off exponentially, 



-xp-x 



p — 



-x(m+p 2 L /2m) 



x (d-2)/2 



(5) 



We are interested not in the familiar exponential fall-off, 
but in power law prefactor of l/x^ 2 ^ 2 . While derived 
in free field theory, it is valid generally. It arises from the 
fluctuations of a point particle, in the directions perpen- 
dicular to the direction along which we measure it. 

In perturbation theory, the Polyakov loop is near unity, 
and can be expanded in powers of Aq . Assuming for the 
purposes of discussion that the Ao field is purely static, 
we obtain 



|^tr(^) 



(6) 



Thus the real part of the Polyakov loop starts out as 
a coupling to two -Ao's, and the imaginary part of the 
Polyakov loop, to three Ao's. Thus to lowest order in 
perturbation theory , 



(t r (x)£ r (0)) - (£)' 



cxp( — 'Zuijjx) 



exp(— 3m,Dx) 



(7) 



(8) 



Here mo is the Debye mass; to lowest order in per- 
turbation theory, mr) — gT . Each result is simply the 
product of either two or three Aq propagators over dis- 
tances x 3> 1/mjj. Thus, if one ignores the difference 
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in prefactors, and defines m r and rtii from the exponen- 
tial falloff of (0) and (@), respectively, we conclude that 
mi/m r = 3/2. These estimates could be altered by mix- 
ing with the purely magnetic sector p] , which we discuss 
at the end of the paper. 

We now turn from the perturbative regime, valid at 
very high temperatures, to temperatures near T c . In- 
stead of dealing with Aq, we construct an effective La- 
grangian directly for the Polyakov loop itself The 
theory must be invariant under the global Z(3) symme- 
try of (^). Implicitly, we assume that the deconfming 
transition for three colors is so weakly first order that 
it is appropriate to classify the effective Lagrangian ac- 
cording to a renormalization group analysis for a second 
order transition. At large distances, x ^> 1/T, the theory 
is effectively three dimensional; by the renormalization 
group, all relevant and marginal operators include terms 
up to sixth order in £. The most general, renormalizable 
Lagrangian in three dimensions, invariant under a global 
Z(3) symmetry, is 



h 
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i 3 + {ty 



e 3 + (i*y 



(9) 



The terms which only involve powers of |^| 2 , proportional 
to ~ 62 1 64 and be , are each each invariant under a global 
symmetry of U(l), £ — > exp(i6)£ for arbitrary 8. The ki- 
netic term is also U(l) invariant. For reasons to become 
clear shortly, we include a factor of wave-function renor- 
malization, Zw , but suppress all other renormalization 
constants for mass and coupling constants. The cubic 
term, ~ 63 £ 3 , and the pentic term, ~ bM\ 2 £ 3 , are only 
invariant under the Z(3) symmetry of (0). Of the two 
hexatic interactions, that ~ be is U(l) invariant, while 
that ~ be£ e is invariant not only under Z(i), but under 
a larger symmetry of Z(6), £ — > exp(iri/3)£. 

It is natural to require that the potential is bounded 
from below for large values of I. At large £, the hexatic 
terms dominate over all other couplings. For real £, this 
gives 



b 6 > 



(10) 



One can also choose £ = exp(z7r/6)£o, with real £q\ then 
the term proportional to be drops out, giving the condi- 
tion 



be > 



(11) 



We perform a mean field analysis, where all coupling 
constants are taken as constant with temperature, ex- 
cept for the mass term, ~ t> 2 2 ■ About the transition, 



condensation of £ is driven by changing the sign of 62 
The sign of 63 is determined by the following. At 
very high temperature, the favored vacuum is perturba- 
tive, with £0 w 1, times Z(3) rotations. We then choose 
63 > so that in the Z(2>) model, there is always one 
vacuum with a real, positive expectation value for £q. 

We also choose 64 > 0. This is not a matter of conven- 
tion. With the above potential, in mean field theory there 
are two ways of obtaining a first order transition. One is 
through a non-zero, cubic coupling. The other is if the 
four point coupling is negative. It is natural, however, 
to assume that 63,65 ^ and 64 > 0. For example, if 
63 = 65 = 0, then the unbroken symmetry of the effective 
Lagrangian is not Z(3), but Z(6). While one can then ob- 
tain a first order transition by choosing 64 to be negative, 
it would also imply a six-fold degeneracy between vacua 
in the broken symmetry phase. Such a degeneracy is not 
observed, even approximately, by Lattice simulations in 
the deconfined phase Consequently, we assume that 
63,65 7^ 0, and that 64 is positive. This also agrees with 
results from the Polyakov Loop Model (^^] . 

While the value of the pentic coupling constant is not 
well fixed, but it cannot dominate the hexatic couplings. 
If the hexatic interactions vanish, 6g = 6g = 0, then 
with £ — exp(i8 /3)£q, the pentic term is negative for 
ir/2 < 9 < 3ir/2, and the potential is unbounded from 
below. Consequently, we consider the pentic interactions 
as a small perturbation on the other couplings. 

This is an effective Lagrangian in three dimensions, 
so that £ has dimensions of mass 1 / 2 ; then the coupling 
constants 6„ have dimensions of mass^ 6 ^™^ 2 . In particu- 
lar, be and 6g are dimensionless. To make the equations 
somewhat simpler, we choose the overall mass scale so 
that 64 = 1; expressions for arbitrary values of 64 can 
be derived from the following by an obvious rescaling. 
We remark that in the Polyakov Loop Model |^-^[, the 
dimensions of £, and of the coupling constants, are all 
made up by powers of temperature. Near T c , however, 
this doesn't matter. 

We now compute the masses for the real and imaginary 
parts of the Polyakov loop. As stated, we choose the 
vacuum expectation value of £, (£) — £q, to be real. The 
equation for £q is just the first derivative of the potential 
for £, 

- h£ - 63 £ 2 + £l + br o £ 4 + {be + be)£l = . (12) 

Computing second derivatives of the potential, the mass 
squared for the real part is: 

m 2 r = - 6 2 - 2b 3 e + S£ 2 + 4b 5 £ 3 + 5(b 6 + b 6 )£ 4 , 

(13) 

while that for the imaginary part is: 



mf = - 6 2 + 26 3 £ + t - 7 b 5 ^ + ih - 26 6 )^ 

5 



(14) 
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The two-point function of Polyakov loops is then, 

2 exp(— m r x) 



(£ r (x)£ r (o)) - (ey 



exp(— mix) 



(15) 



(16) 



Notice that unlike @) and (|), the power law prefactor 
is in each case just ~ 1/x, from the exchange of a single 
^- field. Above T c , rrii ^ m r . 

The global Z(3) symmetry is unbroken below T c : £q = 
0. With our conventions, 62 < 0; the mass squared 
for l T and £i are positive and equal. Thus in our Z('i) 
model, below T c the two-point functions of Polyakov 
loops remain as in ([ll]) and ( |l6[) at large distances, with 
771 i = m r . We remark that the prefactor of ~ l/ir, for 
T < T c and x 3> 1/T, agrees with calculations in Nambu- 
Goto string models |10|. This is true for any number of 
space-time dimensions: the transverse motion of a point 
particle, and the transverse fluctuations of a world-sheet, 
both give the same prefactor of ~ l/x (d ~ 2)/2 . This is 
an important consistency check on models of Polyakov 
loops, such as ||, below T c . 

When 63 or 65 are nonzero, the transition is necessarily 
of first order. Henceforth, we restrict ourselves to the 
critical temperature, T c . This occurs when the nontrivial 
minimum is degenerate with the trivial minimum; i.e., 
when 



h -l£ 2 
7. *n 



J 3 «3 



4 £ ° 



5 to 



= . 
(17) 



This then becomes an equation which determines the 
value of b 2 (T+). 

If we compute each mass individually, m r and m, , then 
there is a complication which we have set aside. The 
kinetic term for £ does not have a standard normalization; 
at tree level, it equals 3/g 2 , where g is the QCD coupling 
constant ||. At one loop order, Wirstam || finds that 
this wave function renormalization is 



9{T) 2 



1 



.08 



g~_ 

' ill 



(18) 



Consequently, the masses measured on the Lattice are ac- 
tually 1/s/Zw of the expressions above. Thus the overall 
scale of m r and rrii depend upon the value of g 2 (T), which 
is not easy to extract directly from Lattice data. 

To avoid the dependence on Zyy, instead we compute 
the ratio of mi/m r , in which Zyy completely drops out. 
It is this ratio which we think would be especially inter- 
esting to compute on the Lattice. 

These equations can be solved numerically for arbi- 
trary values of 65, be, and be- It is illustrative to first 
solve for the case in which 65 = be = be = 0. For arbi- 
trary b 2 , 



to = g I 6 3 + 



46 2 



(19) 



With this £0, at the transition temperature, 62 = ~26g/9, 
and mi/m r = 3. This ratio changes as 65, 6g, and b$ 
increase from zero. To lowest nontrivial order in these 
coupling constants, 



1 



15 



bxb< 



■3 "5 



27 



bl 4& 6 + 5b, 



(20) 



For positive values of b§ and be (remember be is neces- 
sarily positive), the ratio Wi/m r decreases from 3 as 65, 
be, and be increase. It is also easy to obtain mi/m r > 3, 
though. The coupling 65 could be negative; even if 65 = 0, 
when be = —be, etc. 

One can also solve for large values of be and be- Taking 
65 = 0, in this case 




be + be 



(21) 



The combination of be and be on the right hand side is 
necessarily a ratio of positive numbers in all cases. 

Working up from T c , it is also possible to compute 
mi/m r within the context of the Polyakov Loop Model 
(PLM) H . In the PLM, the pressure is determined by the 
potential of £ times T A ; for constant values of 63, 64, 65, 
be, and be, the temperature dependence of the pressure 
then fixes 62 (T) . Consequently, this also determines how 
mi/m r changes with temperature. This appears to be a 
crucial and unambiguous test of the PLM. For example, 
if we take be = be = 0, then at very high temperatures, 
where £ « 1, @ fixes b 2 = 1 - 63, and @ 



rrij 



3b 3 



(22) 



To agree with the perturbative result asT— > 00, = 3/2, 
fixes &3 = 6/7. This is remarkably close to the value 
found from a fit to the pressure near T c , which is 63 w 0.9 
ijljj . A more careful analysis, including the effects of 
pentic and hexatic interactions, is necessary in order to 
make a detailed comparison. 

To summarize, naively in the perturbative regime at 
very high temperature, m,i/m r = 3/2 0. If pentic and 
hexatic terms in £ are a small perturbation on the quar- 
tic terms, then as T — > T c , mi/m r should increase; at 
T c , how close this ratio is to 3 then indicates how im- 
portant these other coupling constants are. Conversely, 
if 77ii/m r does not increase significantly as the tempera- 
ture decreases from 2T C down to T+, then while it might 
be able to obtain a fit with large values of be and be, to 
us this would seem unnatural, and disfavor the present 
analysis. 
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Our discussion neglects mixing between the electric 
sector and the static, magnetic sector J§|. While the for- 
mer might be perturbative, the latter never is. Static, 
magnetic glueballs couple to each other as a three- 
dimensional gauge theory, with coupling constant g 2 T ' . 
Hence the mass of any static, magnetic glueball is some 
number times g 2 T ' . In weak coupling, then, as the Debye 
mass is mu ~ gT, electric fields are much heavier than 
magnetic glueballs. The real part of the Polyakov loop, 
~ tr(Ag), can then mix with spin-zero glueballs, such as 
~ tr(G^ ). Also, the imaginary part of the Polyakov loop, 
~ tr(Ag) can mix with the spin-zero glueball ti(AoGfj). 
(There is also a spin-one glueball, e 13 tr(AoGjk). How- 
ever, the coupling of such a spin-one state to any spin- 
zero state must involve one derivative, and so vanish at 
zero momentum. Spin here refers to a classification in 
three dimensions.) 

Simulations from temperatures of ~ 2T C on up have 
been done by Hart, Laine, and Philipsen |12|; for related 
work, see []l3| . They find that the magnetic sector is not 
lighter than the electric until temperatures T > 10 2 T C . 
Unexpectedly, they do not see significant mixing between 
states such as tr(Ag) and tr(Gy). For example, at ex- 
tremely high temperatures of 10 n T c , the mass of tr(Ag) 
is approximately equal to that of tr(AQ); both masses are 
about ss 1.9 times the mass of the lightest static, mag- 
netic glueball, tr(G?y). The lightest electric state is the 
spin-one glueball, e ljk tr(A Gjk); we note the (presum- 
ably meaningless) numerical coincidence that its mass is 
= 1.5 of the lightest static, magnetic glueball to within 
0.1%, although the quoted statistical errors are w ±6%. 

At temperatures of ~ 2T C , the situation is completely 
reversed: then the electric sector is lighter than that for 
the static, magnetic sector. The lightest state is that for 
the real part of the Polyakov loop, ~ tr(-Ag). Relative 
to this state, the mass of the spin-zero electric state, 
tr(j4g), is « 1.81; that of the spin-one electric state, 
e 13 tr(AoGjk), is « 1.76; that of the lightest static, mag- 
netic glueball, tr(G^), is w 2.4. 

Thus at 2T C , the ratio mi/m r « 1.81. It would be very 
interesting to know how rrii/m r changes at T : 2T C — ► T c \ 
we are unaware of data on mi in this range. Nevertheless, 
at least we can safely neglect mixing of the (light) electric 
sector with the (heavy) magnetic sector. 

At present, most Lattice data is for the two point func- 
tion of complex conjugate Polyakov loops, (£(x)£*(0)} — 
\(£}\ 2 - From either the Z(3) model, or perturbation the- 
ory, in the deconfincd phase it appears that m., is always 
greater than m r . Thus without extracting excited states, 
this correlation function only gives us information on m r , 
and not on to*. 

One can also decide whether an effective model of 
Polyakov loops applies by seeing if the power law pref- 
actor is 1/x, as in jl5|), or 1/x 2 , as in (0). From Lattice 
simulations in 3 + 1 dimensions by Kaczmarek, Karsch, 
Laermann, and Lutgemeier fL4| , the power law prefactor 
is very near unity only near T c , and increases to ~ 1.4 



by 2T C . Thus our effective model does appear to hold in 
at least some limited region near T c . As the exchange of 
any massive state generates a prefactor of ~ 1/x, instead 
of varying the prefactor, to us it seems more natural to 
take a sum over several states, each a prefactor of ~ 1/x. 
We do admit that this is only true in mean field the- 
ory, as (nearly) critical fluctuations introduce nonzero 
anomalous dimensions, rj, which change the prefactor to 
~ 1/x 1- ''. In three dimensions, however, such anoma- 
lous dimensions r) are usually very small; rj < .1 for cither 
Z{2) or 0(2) symmetry groups |l5| ]. 

It is also of interest to consider the transition in 2+1 di- 
mensions. While mean field theory suggests a first order 
transition, fluctuations in two dimensions are so strong 
that the transition is in fact of second order. Then our 
mean field analysis is completely inappropriate; one is at 
the symmetric phase at T c , with = m r . However, cor- 
relation functions of both imaginary and real parts of the 
Polyakov loop have been computed above T c by Bialas, 
Morel, Petersson, and Petrov [jl6f . In all cases, the power 
law prefactor is found to be l/yx, which is what it should 
be for single particle exchange in d = 2 + 1 dimensions. 
Bialas et al. have also computed m,i/m r ; working down 
from high temperatures, this ratio increases to w 2.0, 
then seems to decrease. This might arise from a trade- 
off between mean field effects, as discussed herein, and 
critical fluctuations near T c , special to 2 + 1 dimensions. 

In passing, we also comment that studies with two col- 
ors are also of interest. For two colors, I is real, and so 
there is only m r . One could test, however, to see whether 
the relation between m r and the pressure, as given by the 
PLM, holds up §. With only one mass to measure, the 
value of m r does depend upon Zyp, and so g (T). Also, 
the effects of critical fluctuations must be incorporated. 
On the other hand, the potential for £ only involves two 
coupling constants, 64, and be; 63 = 65 — = 0. Still, 
one basic test is that in 3 + 1 dimensions, the deduced 
value of g 2 (T) increases as T decreases. Lastly, there 
are speculative reasons for thinking that the width of the 
critical region is greater in 2 + 1 dimensions than in 3 + 1 
dimensions 

We look forward to future measurements of these and 
other quantities on the Lattice, near the critical temper- 
ature. 
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